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This first ever double pulsar system consists of two pulsars orbiting the common center of mass in 
a slightly eccentric orbit of only 2.4-hr duration. The pair of pulsars with pulse periods of 22 ms and 
2.8 sec, respectively, conhrms the long-proposed recycling theory for millisecond pulsars and provides 
an exciting opportunity to study the works of pulsar magnetospheres by a very fortunate geometrical 
alignment of the orbit relative to our line-of-sight. In particular, this binary system represents a 
truly unique laboratory for relativistic gravitational physics. This contribution serves as an update 
on the currently obtained results and their consequences for the test of general relativity in the 
strong-held regime. A complete and more up-to-date report of the timing results will be presented 
elsewhere shortly. 


I. INTRODUCTION 

In this year 2005 we celebrate the work of Albert 
Einstein, remembering his enormous contribution to 
our understanding of nature and the Universe. One of 
the best ways for honoring his work is to point out that 
still today, a hundred of years later, hundreds of scien¬ 
tists around the world are deeply involved in searching 
for the limits up to which his centennial theory of gen¬ 
eral relativity (GR) can be applied. To date GR has 
passed all observational tests with flying colours. Still, 
it is the continued aim of many physicists to achieve 
more stringent tests by either increasing the precision 
of the tests or by testing different aspects. Some of 
the most stringent tests are obtained by satellite ex¬ 
periments in the solar system. One must not forget, 
however, that these solar-system experiments are all 
made in the gravitational weak-field regime and that 
they will never be able to provide tests in the strong- 
held limit where deviations from GR may appear more 
clearly or even for the hrst time (see e.g. |l|)- This 
strong-held regime is best explored using radio pul¬ 
sars. 

Pulsars, highly magnetized rotating neutron stars, 
are unique and versatile objects which can be used to 


study an extremely wide range of physical and astro- 
physical problems. Beside testing theories of grav¬ 
ity one can study the Galaxy and the interstellar 
medium, stars, binary systems and their evolution, 
plasma physics and solid state physics under extreme 
conditions. In these proceedings, we will present such 
applications for gravitational physics made possible 
by the hrst ever discovered double pulsar |3,13 ■ We 
will demonstrate that this rare binary system repre¬ 
sents a truly unique laboratory for relativistic gravity. 
We will present an update on the currently obtained 
timing results and their consequences for tests of GR. 
A complete and more up-to-date report of the timing 
results will be presented elsewhere shortly. 


II. THE DOUBLE PULSAR 

Our team discovered the 22.8-ms pulsar 
J0737—3039 in April 2003 Q] in an extension to 
the hugely successful Parkes Multi-beam survey 
13 . It was soon found to be a member of the most 
extreme relativistic binary system ever discovered: 
its short orbital period {Pb = 2.4 hrs) is combined 
with a remarkably high value of periastron advance 
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{(jj = 16.9deg yr“^, i.e. four times larger than for 
PSR B1913+16!) that was measurable after only a 
few days of observations. The system parameters 
predict that the two members of the binary system 
will coalesce on a short time scale of only ~ 85 Myr. 
This boosts the hopes for detecting a merger of two 
neutron stars with first-generation ground-based 
gravitational wave detectors by a factor of several 
compared to previous estimates based on only the 
double neutron stars B1534-I-12 and B1913-I-16 @,S|- 

In October 2003, we detected radio pulses from the 
second neutron star when data sets covering the full 
orbital period were analysed Q . The reason why sig¬ 
nals from the 2.8-s pulsar companion (now called PSR 
J0737—3039B, hereafter “B”) to the millisecond pul¬ 
sar (now called PSR J0737—3039A, hereafter “A”) 
had not been found earlier, became clear when it was 
realized that B was only bright for two short parts of 
the orbit. For the remainder of the orbit, the pulsar B 
is extremely weak and only detectable with the most 
sensitive equipment. The detection of a young com¬ 
panion B around an old millisecond pulsar A confirms 
the evolution scenario proposed for recycled pulsars 
(e.g. i 0) and made this already exciting system 
sensational, providing a truly unique testbed for rela¬ 
tivistic gravity. 

III. STRONG-FIELD TESTS OF GENERAL 
RELATIVITY 

Since neutron stars are very compact massive ob¬ 
jects, the double pulsar (and other double neutron 


star systems) can be considered as almost ideal point 
sources for testing theories of gravity in the strong- 
gravitational-field limit. Tests can be performed when 
a number of relativistic corrections to the Keplerian 
description of an orbit, the so-called “post-Keplerian” 
(PK) parameters, can be measured. For point masses 
with negligible spin contributions, the PK parameters 
in each theory should only be functions of the a priori 
unknown neutron star masses and the well measurable 
Keplerian parameters. 


With the two masses as the only free parameters, 
the measurement of three or more PK parameters 
over-constrains the system, and thereby provides a 
test ground for theories of gravity. In a theory that 
describes a binary system correctly, the PK param¬ 
eters produce theory-dependent lines in a mass-mass 
diagram that all intersect in a single point. 


As A has the faster pulse period (and is bright 
throughout the entire orbit apart from a ~ 27-s eclipse 
at superior conjunction), we can time A much more 
accurately than B and measure precise PK parame¬ 
ters for A’s orbit. In GR, the five most important PK 
parameters are given to first post-Newtonian (IPN, or 
0{v'^ jc^)) order by 0 ; 


" = (g) 

= 7 . 2/3 Mb{Ma + 2 Mb) 

^ ® \ 2 - k ) "" (Ma + Mb)4/3 ’ 

p = 7.5/3 (1 + + je") ^aMb 

^ 5 ® V2^/ (l-e2)7/2 (Ma + Mb)i/3’ 

r = TqMb, 

^ _i /3 (Ma + Mb)^/^ 

~ ® Mb 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


where Pb is the period and e the eccentricity and x 
the semi-major axis measured in light-s of the binary 
orbit. The masses Ma and Mb of A and B, respec¬ 
tively, are expressed in solar masses (Mq). We define 
the constant Tq = GM^jc? = 4.925490947/ts where 
G denotes the Newtonian constant of gravity and c the 
speed of light. The first PK parameter, w, is the eas¬ 


iest to measure and describes the relativistic advance 
of periastron. According to Eqn. ^it provides an im¬ 
mediate measurement of the total mass of the system, 
{M-a + Mb)- The parameter 7 denotes the amplitude 
of delays in arrival times caused by the varying effects 
of the gravitational redshift and time dilation (second 
order Doppler) as the pulsar moves in its elliptical or- 
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bit at varying distances from the companion and with 
varying speeds. The decay of the orbit due to gravi¬ 
tational wave damping is expressed by the change in 
orbital period, Pb- The other two parameters, r and s, 
are related to the Shapiro delay caused by the grav¬ 
itational field of the companion. These parameters 
are only measurable, depending on timing precision, 
if the orbit is seen nearly edge-on. For pulsar A, all 
these quantities have indeed been measured, providing 
a large number of available tests. In fact, in addition 
to tests with these PK parameters, the possibility to 
measure the orbit of both A and B opens up opportu¬ 
nities that go well beyond what is possible with previ¬ 
ously known double neutron stars, as we will describe 
now. 

With a measurement of the projected semi-major 
axes of the orbits of both A and B, we obtain a pre¬ 
cise measurement of the mass ratio, R{Ma, Mb), from 
Kepler’s third law, 

R{Ma, Mb) = Ma/Mb = xbIxa- (6) 

For every realistic theory of gravity, we can expect 
the mass ratio, R, to follow this simple relation [^, 
at least to IPN order. Most importantly, the R 
value is not only theory-independent, but also inde¬ 
pendent of strong-held (self-held) ehects which is not 
the case for PK-parameters. This provides a stringent 
and new constraint for tests of gravitational theories 
as any combination of masses derived from the PK- 
parameters must be consistent with the mass ratio 
derived from Kepler’s 3rd law. With hve PK parame¬ 
ters already available, this additional constraint makes 
the double pulsar the most overdetermined system to 
date where the most relativistic effects can be studied 
in the strong-held limit. 


IV. TIMING OF THE DOUBLE PULSAR 

Our observations already provide measurements for 
all hve PK parameters listed above. This includes 
a measurement of an orbital decay of the binary or¬ 
bit due to gravitational wave emission with a rate of 
7mm/day. As indicated earlier, we can use these re¬ 
sults to test GR in a very elegant way Q . The unique 
relationship between the two masses of the system pre¬ 
dicted by GR (or any other theory) for each PK pa¬ 
rameter can be plotted in a diagram showing the mass 
of A on one axis and B on the other. We expect all 
curves, including that of the mass ratio R, to intersect 
in a single point if the chosen theory (here GR) is a 
valid description of the nature of this system. Such 
tests have been possible to date in PSR B1913-I-16 
(e.g. M) and for PSR B1534-fl2 (e.g. M)- How¬ 
ever, in neither of these systems were so many curves 
available as for the double pulsar system for which we 
derive a Ma — Mb plot as shown in Fig.^ 


It turns out that, as another stroke of luck, we 
are observing the system almost completely edge-on 
which allows us to determine a Shapiro delay to very 
high precision. It also means, however, that we can 
probe a pulsar magnetosphere for the very first time 
using a background beacon. Results of the mag- 
netospheric interactions between A and B and 
the eclipse of A’s signal at its superior conjunction 
13, [ 13 , nil have been presented elsewhere and must 
be considered in search for a possible contamination 
of the timing data and hence a violation of our as¬ 
sumption that we deal with a “clean” system of point 
sources. 

While all studies so far confirm the cleanness of the 
system, we also have to consider that the times-of- 
arrival (TOAs) are obtained with a standard “tem¬ 
plate matching” procedure that involves a cross¬ 
correlation of the observed pulse profile with high 
signal-to-noise ratio template (e.g. [13 )• Any change 
in the pulse profile could therefore lead to systematic 
variations in the measured TOAs. For this reason it 
was necessary to perform detailed studies of the pro¬ 
files of A and B and to investigate any possible profile 
evolution with time. Indeed, profile changes on secu¬ 
lar time scales are expected for the double pulsar for 
the following reason. 

In GR, the proper reference frame of a freely falling 
object suffers a precession with respect to a distant ob¬ 
server, called geodetic precession. In a binary pulsar 
system this geodetic precession leads to a relativistic 
spin-orbit coupling, analogous of spin-orbit coupling 
in atomic physics As a consequence, the pul¬ 

sar spins precess about the total angular momentum, 
changing the relative orientation of the pulsars to one 
another and toward Earth. Since the orbital angular 
momentum is much larger than the pulsars’ angular 
momenta, the total angular momentum is practically 
represented by the orbital angular momentum. The 
precession rate depends on the period and the ec¬ 
centricity of the orbit as well as the masses of A and 
B. With the orbital parameters of the double pulsar, 
GR predicts precession periods of only 75 yr for A and 
71 yr for B. 

Geodetic precession should have a direct effect on 
the timing as it causes the polar angles of the spins 
and hence the effects of aberration to change with time 
[^ . These changes modify the observed orbital param¬ 
eters, like projected semi-major axis and eccentricity, 
which differ from the intrinsic values by an aberration 
dependent term, pot entially allowing us to infer the 
system geometry [T^. Other consequences of geode¬ 
tic precession can be expected to be detected much 
sooner and are directly relevant for the timing of A 
and B. These arise from variations in the pulse shape 
due to changing cuts through the emission beam as the 
pulsar spin axes precess. Moreover, geodetic preces¬ 
sion also leads to a change in the relative alignment of 
the pulsar magnetospheres, so that the visibility pat- 
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TABLE I: Observed and derived parameters of PSRs J0737—3039A and B. Standard errors are given in parentheses after 
the values and are in units of the least significant digit(s). 


Pulsar 

PSR J0737-3039A PSR J0737-3039B 

Pulse period P (ms) 

22.699378556138(2) 2773.4607474(4) 

Period derivative P 

1.7596(2) X 10-18 0.88(13) x lO"!® 

Epoch of period (MJD) 

52870.0 

Right ascension q (J2000) 

07® 37“51=.24795(2) 

Declination S (J2000) 

-30°39'40".7247(6) 

Orbital period (day) 

0.1022515628(2) 

Eccentricity e 

0.087778(2) 

Epoch of periastron To (MJD) 

52870.0120588(3) 

Advance of periastron u) (deg yr~^) 

16.900(2) 

Longitude of periastron lj (deg) 

73.805(1) 73.805 -1- 180.0 

Projected semi-major axis x = asini/c (sec) 

1.415032(2) 1.513(4) 

Gravitational redshift parameter 7 (ms) 

0.39(2) 

Shapiro delay parameter s = sin 2 

0.9995(4) 

Shapiro delay parameter r (ps) 

6 . 2 ( 6 ) 

Orbital decay Pf, (10“'^^) 

- 1 . 20 ( 8 ) 

Mass ratio R = 

1.071(1) 


tern or even the profile of B may vary due to these 
changes as well. 

The possibility to observe such phenomena is ex¬ 
citing but requires a careful analysis to exclude any 
impact onto the timing of the two pulsars. Indeed, 
the necessity to check for these effects has delayed the 
publication of the currently final timing results of the 
double pulsar. The analysis is now almost complete 
and results will be published shortly. A preliminary 
update is given in Table 1 and discussed below. 


V. PRESENT RESULTS 

The study of the profile evolution of A did not 
lead to the detection of any profile change over a pe¬ 
riod of 15 months. This present non-detection greatly 
simplifies the timing of A but does not exclude the 
possibility that changes may not happen in the fu¬ 
ture. While the effects of geodetic precession could 
be small due to a near alignment of pulsar A’s spin 
and the orbital momentum vector, the results could 
also be explained by observing the system at a par¬ 
ticular precession phase. While this case appears to 
be relatively unlikely, it must not be excluded as such 
a situation had indeed occurred for PSR B1913-I-16 
[^. Indeed, a modelling of the results suggests that 
this present non-detection of profile changes is con¬ 
sistent with a rather wide range of possible system 
geometries. One conclusion that can be drawn, how¬ 
ever, is that the observations are inconsistent with the 
large profile changes that had been predicted by some 
models m. 

In contrast to the results for A, similar studies of 


the profile evolution of B reveal a clear evolu¬ 
tion of B’s emission on orbital and secular time-scales. 
The profile of B is changing with time, while also the 
light-curves of B (i.e. the visibility of B versus orbital 
phase) undergo clear changes. These phenomena may 
be caused by a changing magnetospheric interaction 
due to geometry variations resulting from geodetic 
precession. In any case, these changes require sophis¬ 
ticated timing analysis techniques and the preliminary 
results obtained with this techniques are listed in Ta¬ 
ble 1, while final results will be published shortly. 

The present timing results already indicate that the 
proper motion of this system is surprisingly small. 
While a significant measurement of a proper motion 
vector via pulsar timing will be available shortly, the 
present limit suggests a systemic velocity of less than 
30 km/s for a dispersion measure distance of 600 pc 
. [33 While such a small velocity may be indicative 
of a small kick imparted onto B during its supernova 
explosion |^, other studies find this limit still to be 
consistent with a kick of average magnitude |26|. In 
any case, such a small velocity is good news for tests of 
GR. Usually, the observed value of Pb is positively bi¬ 
ased by an effect known as secular acceleration arising 
from a relative motion and acceleration of the system 
(e.g. 113 )■ Computing the magnitude of this observa¬ 
tional bias using the obtained limit on the proper mo¬ 
tion, however, suggests that the contribution is much 
less than 1%, so that the orbital decay measurement 
will be available for another precise GR test. 

Finally, scintillation measurements have recently 
suggested an orbital inclination angle that is ex¬ 
tremely close to 90 deg (i.e. within (0.29 ± 0.14) deg) 
This measurement appears to be inconsistent 
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Mass A (Mg^jj) 

FIG. 1: ‘Mass-mass’ diagram showing the observational constraints on the masses of the neutron stars in the double 
pulsar system J0737-3039. The shaded regions are those that are excluded by the Keplerian mass functions of the 
two pulsars. Further constraints are shown as pairs of lines enclosing permitted regions as given by the observed mass 
ratio and PK parameters shown here as predicted by general relativity. Inset is an enlarged view of the small square 
encompassing the intersection of these constraints (see text). 


with results from the timing observations and the 
measurement of the Shapiro delay parameter, s = 
sini, which suggest an inclination angle that is close 
but significantly different from 90 deg. One should 
note that the scintillation results are based on corre¬ 
lating the scintillation properties of A and B over the 
short time-span of the orbital motion when they are 
in conjunction to the observer. In contrast, the mea¬ 
surement of the inclination angle from timing mea¬ 
surements results from detecting significant harmonic 
structure in the post-fit residuals after parts of the 
Shapiro delay effect are absorbed in the fit for the 
Romer delay, i.e. the light travel time across the or¬ 
bit. As shown in Figure[21 these structures are present 
throughout the whole orbit, so that the results from 
timing measurements may be expected to be more re¬ 
liable. We are currently studying the origin of this 
apparent inconsistency between these two methods, 
checking both any contamination of the Shapiro de¬ 
lay measurements and effects influencing the scintilla¬ 
tion results. An exciting possibility could be that the 


emission of A suffers measurable refraction while prop¬ 
agating through the magnetosphere of B. If that were 
indeed the case, we would have a direct handle onto 
the magneto-ionic properties of B’s magnetosphere for 
the first time. 


Inspecting the results shown in Table I, we can 
take the most precise parameters (i.e. the mass ra¬ 
tio R, the advance of periastron to and the Shapiro 
delay parameter s) to test theories of gravity. As¬ 
suming that GR is the correct theory of gravitation, 
we use Eqn.nto derive the total mass of the system 
and combine it with the observed mass ratio to obtain 
Ma = 1.338 ± O.OOlMo and Mb = 1.249 ± O.OOlMg. 
Using these precisely determined masses we compute 
the Shapiro delay parameter s as predicted by GR 
and compare it to the observed value. We find that 
gGR/gobs _ i.0002lg;QQjg. Hence, GR passes this test 
at the 0.1% level. This is the most stringent test of 
GR in the strong-held limit so far. 
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FIG. 2: The effect of the Shapiro delay caused by the 
gravitational potential of B seen in the timing residuals 
of A. The timing residuals obtained by fitting all model 
parameters shown in Table 1 except the Shapiro delay pa¬ 
rameters r and s. The left-over structure represents the 
higher harmonics of the Shapiro delay that are unabsorbed 
by fits to the Keplerian parameters. 

VI. FUTURE 

In the near and far future, the precision of the de¬ 
termined parameters will increase further, simply by 
the available longer time span but also by the po¬ 
tential employment of better instrumentation. In a 
few years, we should be able to measure additional 
PK parameters, including those which arise from a 
relativistic deformation of the pulsar orbit (resulting 
in angular and radial orbital eccentricities) and those 
which find their origin in aberration effects and their 
interplay with geodetic precession (see i). On sec¬ 
ular time scales we will even achieve a precision that 
will require us to consider post-Newtonian terms that 
go beyond the currently used description of the PK 
parameters. Indeed, the equations for the PK param¬ 
eters given earlier are only correct to lowest PN order. 
However, higher-order corrections are expected to be¬ 
come important if timing precision is sufficiently high. 
While this has not been the case in the past, the dou¬ 
ble pulsar system mw allow measurements of these 
effects in the future Q . 

One such effect involves the prediction by GR that, 
in contrast to Newtonian physics, the neutron stars’ 
spins affect their orbital motion via spin-orbit cou¬ 
pling. This effect would be visible most clearly as a 
contribution to the observed ti; in a secular |l 7l| and 
periodic fashion [^. For the J0737—3039 system, the 
expected contribution is about an order of magnitude 
larger than for PSR BI9I3-I-16, i.e. 2 x 10“"^ deg yr“^ 
(for A, assuming a geometry as determined for PSR 
BI9I3-I-16 [13). As the exact value depends on the 
pulsars’ moment of inertia, a potential measurement 
of this effect allows the moment of inertia of a neutron 


star to be determined for the first time [^. 

If two parameters, e.g. the Shapiro parameter s and 
the mass ratio R, can be measured sufficiently accu¬ 
rate, an expected Wexp can be computed from the in¬ 
tersection point. This value can be compared to the 
observed value Wobs which is given by (see [^) 

Wobs = WlPN [l + Ad;2PN ~ ~ 9^^^So\ 

(7) 

where the last two terms represent contributions from 
the pulsar spin. In these terms, are geometry 

dependent factors whilst Auj^^ arise from relativistic 
spin-orbit coupling, formally at the IPN level. How¬ 
ever, it turns out that for binary pulsars these effects 
have a magnitude equivalent to 2PN effects |3, so 
that they only need to be considered if oj is to be 
studied at this hi gher level of approximation. We find 
Awso cc I/PM"^ 1^, so that with precise masses M 
the moment of inertia I can be measured and the neu¬ 
tron star “equation-of-state” and our understanding 
of matter at extreme pressure and densities can be 
tested. 

The dependence of AAso on the spin period P sug¬ 
gests that only a measurement for pulsar A can be 
obtained. It also requires that at least two other pa¬ 
rameters can be measured to a similar accuracy as to. 
Despite being a tough challenge, e.g. due to the ex¬ 
pected profile variation caused by geodetic precession, 
the prospects are promising. Simulations indicate that 
a few years of high precision timing are sufhcient. 


VII. SUMMARY & CONCLUSIONS 

With the measurement of five PK parameters and 
the unique information about the mass ratio, the PSR 
J0737—3039 system provides a truly unique test-bed 
for relativistic theories of gravity. So far, GR also 
passes this test with flying colors. The precision of 
this test and the nature of the resulting constraints go 
beyond what has been possible with other systems in 
the past. The test achieved so far is, however, only the 
beginning of a study of relativistic phenomena that 
can be investigated in great detail in this wonderful 
cosmic laboratory. 
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